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Abstract 

The SU(3) extension of the linear sigma model is employed to elucidate the 
effect of including strangeness on the formation of disoriented chiral conden- 
sates. By means of a Hartree factorization, approximate dispersion relations 
for the 18 scalar and pseudoscalar meson species are derived and their self- 
consistent solution makes it possible to trace out the thermal path of the two 
order parameters as well as delineate the region of instability within which 
spontaneous pair creation becomes possible. The results depend significantly 
on the employed sigma mass, with the highest values yielding the largest re- 
gions of instability. An approximate solution of the equations of motion for 
the order parameter in scenarios emulating uniform scaling expansions show 
that even with a rapid quench only the pionic modes grow unstable. Never- 
theless, the rapid and oscillatory relaxation of the order parameters leads to 
enhanced production of both pions and (to a lesser degree) kaons. 
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I. INTRODUCTION 



High-energy collisions of atomic nuclei may produce extended regions of space within 
which chiral symmetry is temporarily nearly restored. The subsequent non-equilibrium 
relaxation towards the normal vacuum may then produce disoriented chiral condensates 
(DCC), a coherent oscillation of the pion field that is manifested through anomalous pion 
multiplicity distributions and associated enhancements of electromagnetic processes 

(for reviews of the topic, see Refs. [p|-pT|). In order to assess the prospects for such 
DCC phenomena to actually occur and be detectable above the background, it is necessary 
to perform extensive dynamical calculations. Most dynamical studies have been carried out 
within the framework of the SU(2) linear a model which incorporates the a and ir degrees 
of freedom by means of an 0(4) field 

While the omission of the strangeness degree of freedom is justified in many nuclear 
studies, since the mass of the strange quark is large compared with the temperatures typically 
encountered, it may not be appropriate in the DCC context because of the high temperature 
required for the approximate restoration of chiral symmetry. We therefore find it of interest 
to extend the treatment to the strange sector and assess the effect of the DCC physics. 

The resulting extended theory has a richer structure, as two order parameters are now 
involved, the light quark condensate and the strange quark condensate, and the group of 
mesonic excitations has expanded from SU(2) to SU(3). Moreover, whereas the usual SU(2) 
treatment considers only the four fields representing a and 7r excitations, the SU(3) model 
includes the eighteen meson fields of the scalar and pseudoscalar nonets. 

The incorporation of strangeness causes significant changes, as will be discussed for 
thermal equilibrium as well as the coupled dynamics of the two order parameters. Moreover, 
we explore the possibility that any of the mesonic quasi-particle modes become unstable 
during the evolution. The study is carried out within the framework of the linear o model. 
Although its SU(3) version has been known for a long time, to our knowledge it has not yet 
been utilized in connection with disoriented chiral condensates. 

The presentation is organized as follows. In Sect. |I| we introduce the SU(3) linear a model 
and fix the parameters to the meson masses in the mean-field approximation. The regions 
of instability of the various meson species are then delineated in the plane of the two order 
parameters and at zero temperature (Sect. pTT[ ) . Subsequently, in Sect. [TV], we study how the 
properties develop with temperature. After having elucidated the equilibrium features, we 
then, in section [V], consider various idealized dynamical scenarios that emulate what may be 
experienced in the rapidly expanding collision zone following a high-energy nuclear collision. 
Finally, we briefly summarize our results. 



II. SU(3) LINEAR a MODEL 



The linear a model was first introduced by Gell-Mann and Levy p7 |. In a subsequent 



work, Levy extended this model from SU(2) to SU(3) |28| based on the success of the 
Eightfold Way |29| (for a review of the early work done in SU(3) chiral models, see Ref. IfjCfl ). 
It turns out that the model successfully accounts for the masses of the pseudoscalar and scalar 
nonet, especially for the 77—77' splitting by including a term that breaks U(3) xU(3) symmetry 
but conserves the SU(3) xSU(3) symmetry (see, for example, Refs. p8|j3l| -p4]| ) . The data for 
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pion-pion and pion-kaon scattering can also be described |35|j36[| and the tree-level predictions 
are in agreement with Weinberg's predictions j37j. Nevertheless, the nonlinear realization 
was preferred over the linear one during the subsequent decades. 

The linear description of chiral symmetry by means of the scalar mesons as explicit 
degrees of freedom is supported by the recent confirmation of the a meson [55||, which has 
led to the reappearance of the a meson in the compilation by the Particle Data Group [|39 



albeit with a very uncertain mass in the range of 400 — 1200 MeV. Interestingly, the recent 
fit to the scalar mesons by Tornqvist |^0[ can in fact be interpreted as a model unitarization 
based upon U(3)xU(3) chiral symmetry ET . 



The SU(3) linear a model was formulated in the original paper by Levy with the baryon 

HI!. The 



octet |p8|| . Vector mesons can be included also as massive Yang-Mills fields 
trace anomaly can be effectively accounted for by a glueball field ||43|| . The glueball field 
might be relevant for the chiral dynamics as pointed out most recently for the SU(2) linear 
a model with quarks [31|. We leave the study of the effects of the glueball field on the DCC 
formation in SU(3) for future study. This model, including baryons, vector mesons, and a 
dilaton field, has recently been applied to describe baryon masses as well as nuclear matter 
and hypernuclei pESfl . 

In the following we discuss the model for pseudoscalar and scalar mesons only, as we are 
interested in meson-dominated environments, such as those expected to be produced in the 
central rapidity region by relativistic nuclear collisions at RHIC. Due to their higher masses, 
vector mesons will become important only at higher temperatures and will therefore not be 
considered here. 

The model used has a strong connection to QCD. The order of the chiral phase tran- 
sition of QCD was discussed using an effective Lagrangian for arbitrary number of flavors 
by Pisarski and Wilczek K| using the linear SU(N) a model which is the most general 
renormalizable Lagrangian consistent with the symmetries of QCD. The chiral critical point 
was examined in [17] using the same effective chiral Lagrangian for 2+1 flavors. We will 
employ the same effective Lagrangian in the present study. 



A. The SU(3) meson states 

For the SU(2) version of the model, there are only four mesonic modes, the a isosinglet 
and the 7r isotriplet, whereas the expanded SU(3) chiral model contains 18 different mesons, 
the nine pseudoscalar and the nine scalar mesons of the respective nonets. The building 
block of the chiral SU(3) symmetric Lagrangian is then the corresponding matrix M, 



M 



E + ill = -jEY.K{<y a + i^a) 

V 1 a=0 



(1) 



Here a a and 7i a denote the U(3) components of the scalar and pseudoscalar nonets, respec- 
tively, and A a denotes the standard Gell-Mann matrices (including Ao = \f\l)- 
The explicit representation in terms of the physical states is 

( 75^° + 75^ + 7^° 

71 

K~ 



7T 



n 



V 



-7^° + 75^ + 73^0 
K° 



K + 
K° 



\ 



2) 



T^ 8 + 7^° / 
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where rjo and T]s are mixed states of the 77 and rj' mesons and we note in particular the 
expressions for the pions and kaons in terms of the Gell-Mann matrices, 



/ K~ 



-j= (Ai7Ta ± zA 2 7r 2 ) 

-- A= (A 4 7T4 ± «A 5 7r 5 ) 



7T 



For the scalar nonet we have 



V 



1 

V2 



a + a 
a 

K *- 



K° 



a 



K 



A 3 7r 3 



-j= (A 6 7T 6 ± iA 7 7T 7 ) 



K *+ 

K*° 

c / 



(3) 
(4) 

(5) 



where we have introduced the two order parameters, a and (, corresponding to the light 
quark and the strange quark condensate, respectively. The ( is in this notation a pure ss 
state, i.e. ideal mixing is assumed. Nevertheless, as we will see, the a and £ can still mix. 



B. SU(3)xSU(3) chiral Lagrangian 



We start from the U(3)xU(3) chirally symmetric renormalizable Lagrangian for scalar 
and pseudoscalar mesons, 



C 



U(3) 



-Tr 



MM f ] - ATr |~MM t MM t J - A' (Tr [MM* 



(6) 



Contrary to the SU(2) case, there are now two coupling terms of fourth order in the fields, 
with the strength parameters A and A'. 

In order to obtain a correct description of the meson masses, it is necessary to include a 
term that breaks the U(3)xU(3) symmetry while still maintaining SU(3)xSU(3) symmetry, 

M + M+ 



C 



SU(3) 



c det Tr 



(7) 



This term splits it and v( which are degenerate in U(3) xU(3), even in the presence of explicit 
symmetry breaking [Q. It is related to the Ua(1) anomaly of QCD and instanton-induced 
interactions (see Ref. and references therein). There are no additional independent terms 
that are invariant under chiral SU(3) xSU(3) transformations, so any higher-order forms can 
be reexpressed as functions of the three interaction terms already introduced. 

The symmetry is explicitly broken by an additional linear term in the Lagrangian, 



5C = f cr + f 8 a 8 = m\j v a + V2m 2 K f, 



K 



V2 



m 2 J, C 



(8) 



which makes the pseudoscalar mesons massive. The second relation follows from the PCAC 
relations for pions and kaons [45,3Cj. The vacuum expectation values for the a and ( fields 
are then fixed by the pion and kaon decay constants, 

:/ 7r = 92.4MeV , (9) 

: V2f K ~ -j^A = 94.5 MeV , (10) 

using the experimental values for the charged meson decay amplitudes, f n ± = 92.4 ± 0.3 
MeV and f K ± = 113.0 ± 1.3 MeV. We note that <r vac and Cvac are remarkably similar. 



a v 



Cv 
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C. SU(2) case 



In the limit of pure SU(2), we consider only the a and 7r fields in the above expressions 
for the meson matrices (0-0). As the ( field is assumed to be zero, the SU(3) term (^) 
vanishes. Moreover, the interaction terms proportional to A and A' are identical because the 
relation 



Tr 



MM t MM t | = \ (Tr [MM*!) (11) 



holds for traceless matrices due to the Burgoyne's identity The matrices are not traceless 
in the general U(3)xU(3) case due to the presence of the singlet states. In SU(2) the 
Lagrangian then reduces to 

£su ( 3Hsu(2) = Adn + ^ 2 (V + vr 2 ) - i(A + 2 A') (a 2 + vr 2 )* + ^m 2 a , (12) 
which is recognized as the SU(2) linear a model, 

£su(2) = ^kin - iA SU (2) (a -2 + Tr 2 - w 2 ) + Umla , (13) 



if one assigns Agu(2) = 2A + 4A' and v 2 = yU 2 /Asu(2) and ignores the constant term |Asu(2)^ 4 . 
However, since the SU(3) case has additional contributions from the strange quark conden- 
sate, the effective parameter values are shifted and hence the values employed in the SU(2) 
linear a model cannot be used as a guideline for SU(3). Consequently, the parameters of 
the SU(3) model will be fixed to the meson masses of the nonets. 



D. Meson masses 

The remaining four parameters of the model, /i 2 , A, A', and c, are fixed by the meson 
masses in vacuum which are given by the eigenvalues of the curvature tensor, 

d<f>id<pj 

evaluated at the equilibrium point. In the mean-field approximation, the only mixings are 
between a and ( and between r\ and rf . The values of the model parameters are constrained 
by the fact that the existence of a stable equilibrium requires M. to be positive definite, i.e. 
the eigenvalues (the physical unmixed masses squared) must be positive. The expressions 
for the masses and the mixing angles are listed in the Appendix. 

The easiest way is to fix the coupling constants to the pion and kaon masses and the 
mass average |(^ 2 + m 2 ,) P^j . The rj-r]' mass splitting is then predicted. Then /i 2 and A' 
are related but can only be fixed by a choice for the scalar meson mass, as discussed in detail 
in Ref. [0. The parameter A' must be larger than —2.85 to ensure a positive value of m 2 . 
The mass parameter /i 2 is only positive for m a > 724 MeV. As already noted, m a is not 



well known, with current estimates ranging between 400 and 1200 MeV Therefore, we 
consider several cases within this range, m a = 400, 600, 800, and 1000 MeV. The parameters 
A, c, and fi are analytically related to the pseudoscalar masses, 
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A 



(ml - \(m 2 v + m 2 v ,)) 


(2Cvac — V^vac^ 


+ 3C 


me (m 2 K ~ ml) 


Kac + 4C4c) (2Cvac - V2a^) 





+ 2ACv 



V^CTvac - 2( v 

/X 2 = 4A' (cr 2 ac + (L) -< + 2ACvac " 2cCvac , (15) 

and A' is then determined by m a . The (unphysical) case cr vac = v^2~Cvac yields degenerate 
masses for strange and non-strange mesons. These parameters are listed in Table | together 
with the masses of the pseudoscalar and scalar meson nonets. 

The pseudoscalar mixing angle turns out to be 8 P = —5°. For comparison, the Gell- 
Mann-Okubo mass formula gives 8 P = —10° while experiment seems to favor 8 P = —20° |39 



Nevertheless, the masses of rj and rj are reproduced at the level of 1-2%. The parameters 
can be changed to get a pseudoscalar mixing angle of — 10?4 for a higher kaon mass and a 
lower kaon decay constant [|47|. A more interesting approach to get a correct mixing angle 



has been recently pointed out in [^9[ where the mass relations were modified by the inclusion 
of the pseudoscalar decay constants. 

The anomaly term (eq. (|7|)) is important for obtaining the correct meson mass spectra, 
especially for giving the r\ its large mass and splitting the tt and r\ masses. In addition, 
the term mixes the 7] and rj' meson states away from the ideal mixing states rj ns and rj s , the 
former being a pure light-quark system and the latter a pure ss state. There are arguments 
that the SU(3) chiral phase transition and the restoration of the Ua(1) chiral symmetry 
does not necessarily coincide (scenario 1 considered by Shuryak [|50|]). Then, a nonvanishing 
anomaly term would be responsible for giving the (tt, a) and (f] ns , a ) meson pairs different 
masses, even in the SU(3) chiral symmetric phase. Their mass splitting is proportional to 
the ( field, hence is related to the finite strange quark mass. Also the r/ ns gets a contribution 
from the anomaly term due to the finite strange quark mass, even in the SU(3) chiral limit, 
but the rj s does not. This could make the r] s lighter than the r] ns , as pointed out by Schafer 
51 1 . On the other hand, the K and K* mesons are degenerate for a vanishing a field, so 
their mass splitting from the anomaly term is related only to the light quark masses. 

We recall that the situation for the scalar nonet is still somewhat controversial (see 
Ref. [^, p. 355] for a discussion of scalar mesons and their assignment to the scalar nonet). 
Weinstein and Isgur suggested, that the / (980) as well as the a (980) are bound KK 
molecules ]52| . The corresponding meson states for the scalar nonet are higher in mass and 
are identified with / (1370) and a (1450). In this picture, the mass of the ( can be chosen to 
be close to / (1370) but then the a is too light (see Table |). On the other hand, Tornqvist 
argues that all the lowest-lying scalar mesons are describable in a quark model and can 
be assigned to the scalar nonet; the more massive scalar states are then excited states. In 
Table [I] we choose the latter picture for comparison. The mass of a (980) is close to its 
experimental value. The ( is much heavier than the /o(980), especially for the case m a = 1 
GeV. The mass of the K* seems to be far off compared to the i^*(1430). Nevertheless, it 
was recently argued that there is evidence for a light scalar k(900) 
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E. Idealized DCC's 



In the SU(2) case, isospin invariance implies that the pion field vector (711,772,773) is 
equally likely to point in any particular direction, tt = (cos (p n sin^^, sin 92^ siniV, cos^). 
If the iso-orientation of the pion field is the same throughout the entire source, then all 
the radiated pions emerge with the same isospin alignment in each separate event. It then 
follows that the resulting neutral pion fraction R n = N n o/N n = cos 2 ^ has the distribution 
P 7T (R n ) = \j2\[R~ K This distribution differs significantly from the narrow peak around 

one third that would be expected if each individual pion had a random isospin orientation 
and it has therefore been proposed as a suitable DCC signal. 

This situation remains the same in the SU(3) case with regard to the three-dimensional 
pion subspace. But, in addition, there is also invariance in the four-dimensional subspace of 
the kaon field (K + , K~, K°, K°), i.e. the direction of the kaon field vector (7r 4 , 7r 5 , 7r 6 , 7r 7 ) 



K = (cos 



ipx sin $k sin ip K , sin ip K sin -d K sin %fj K , cos fix sin ip K , cos ipx ) , (16) 

is equally likely to point in any direction. Thus, by analogy with the pions, it is natural to 
consider the neutral kaon fraction, 

D - N K + N K . 2 . 2 

RK = N K+ + N K - + N K0 + N R0 =1 ~ Sm ^ • 

Under the similar assumption of full alignment of the kaon field in each distinct event, it is 
elementary to calculate the corresponding distribution, 

j rln 

Pk(Rk) = tt^ d?p K sin 2 ip K dd K $md K l d(p K S(R K - 1 + sin 2 d K sin 2 ip K ) (18) 
2n z Jo Jo Jo 

I 1 

= — / dip k sin ipK 6 ] {Rk — cos 2 ipx) [Rk — cos 2 ^x]" 5 (19) 

TT Jo 

= -/ dcosipK [R K - cos 2 i/jr]' 1 = 1- (20) 
n Jo 

Hence, all values of the neutral kaon fraction Rk are equally likely, whereas a statistical 
emission would yield a narrow distribution peaked at one half. As was the case for P n (R 7T ), 
the result for Pk(Rr) is independent of any explicit chiral symmetry breaking, as it relies 
only on the invariance with regard to transformations among the four kaon states. 

Contrary to the neutral pions, neutral kaons may decay by weak interactions into 
hadrons, Kg — > 7r + 7r~, which can be detected in heavy- ion experiments by standard tech- 
niques. Although event-by-event reconstruction is impractical, fluctuation analysis of the 
abundance of neutral and charged kaons has been proven feasible with the present heavy-ion 
experiments at the SPS. An enhancement of the ratio near zero or one, or at least a broaden- 
ing of Pk(Rr), might be visible and would suggest the occurrence of non-equilibrium effects 
of the DCC type. 



III. REGIONS OF STABILITY AT ZERO TEMPERATURE 

In this section we examine the regions of instability for the pseudoscalar mesons in 
the plane of the two order parameters a and (, at zero temperature. This information is 
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particularly relevant in the idealized quench scenario, in which the initially hot system, in 
which chiral symmetry is at least partially restored, is suddenly cooled down. Generally, 
the effective meson masses depend on both condensates and it is thus more complicated to 
map out the occurrence of instabilities. 

Lattice gauge calculations indicate that although the light-quark condensate does not 
entirely vanish in the chirally restored phase it is considerably reduced. The strange-quark 
condensate follows this trend but its reduction is less pronounced due to the higher mass of 
the strange quark |54]]. The different behavior of the strange-quark condensate is ensured 
in the linear SU(3) a model with the inclusion of the ( field. The mean-field expressions 
for the effective meson masses as functions of the order parameters and the thermal field 
fluctuations are given in the Appendix. In the following, we study only the pseudoscalar 
mesons as they are directly measurable in heavy-ion experiments. 



In the mean-field approximation |5q| , the quasi-particles satisfy the Klein-Gordon dis- 
persion relation, io\ = k 2 + m 2 . Thus, whenever a negative value of m 2 occurs, the frequency 

is imaginary for sufficiently small values of the wave number k, leading to an exponential 
growth of the corresponding field component. This phenomenon is familiar from the discus- 
sion of supercritical fields and it manifests itself as an enhancement of the number of soft 
mesons of the affected kind. Therefore it is of interest to delineate the regions of stability 
for the various meson species. 

Figure [l] shows the critical boundaries for the pseudoscalar mesons for the four selected 
values of m a . We see that m 2 is negative in the upper-left part while m 2 is negative in 
the lower-right part. For the relatively small mass of 400 MeV the ir and 77 are the only 
mesons exhibiting instabilities and the corresponding supercritical regions are fairly small. 
In particular, there is no combination of the two order parameters a and ( that would cause 
the kaons to become unstable. These general features remain true also at m a = 600 MeV, 
although the regions of tt and r] instability have expanded considerably. 

This expansion continues steadily as m a is increased further and by 800 MeV the situation 
has changed dramatically, as the model parameter /i 2 has turned positive. Instabilities are 
now possible for all four pseudoscalar mesons and occur simultaneously in a small region 
near the origin where full chiral symmetry is approximately restored. However, it should be 
noted that the instabilities of K, 77, or rj' occur only for considerably reduced values of the 
strange order parameter (. Thus, a reduction of only the light-quark condensate a could 
produce only the standard pionic instability. Moreover, the critical boundaries for the r] and 
K mesons are quite close to one another. 

Finally, for m a = 1 GeV, an imaginary pion mass is reached already after only a slight 
reduction of either o or (. Furthermore, the kaon is now unstable over a wide order- 
parameter region and the critical boundary for rj is again close to the kaon boundary. We 
also note that the critical boundary for rf has expanded considerably so that this meson 
may exhibit instability as well even if chiral symmetry has been only partly restored. 



IV. FINITE TEMPERATURE 



We now study how the properties of the system develops as the temperature is raised. 
Statistical equilibrium of the SU(3) linear a model can be treated in the mean-field approx- 
imation along the lines developed for the SU(2) case [Q. The local field strength is then 
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decomposed into the (smooth) order parameter and the residual fluctuations, 

a — (a) + 5a and C — (0 + • (21) 

The insertion of this decomposition into the general equation of motion then leads to separate 
equations of motion for the order parameters and the individual quasi-particle fields. The 
trilinear fluctuation terms are included by the Hartree factorization technique, 

S^S^jd^k — > (S&Scf) j)5(f) k + (5(f) j 5(f) k }5(f)i + (SfaS^S^j (22) 

where (•) denotes the thermal average, 

«*m - ^ h^i ^lm-i = ^ g k«> (T) • < 23 > 

in the diagonal representation. The sum is over all finite wave numbers k and the associated 



energy is uo k = y mf + k 2 where m; is the appropriate temperature-dependent effective mass 
for the meson species i, as obtained by solving the coupled equations selfconsistently. We 
need to take into account only the thermal fluctuations of ir, K, a, and rj mesons, as they 
are the lightest and thus have the largest field fluctuations. The other mesons (r/, ao, K*, 
() have masses near or above one GeV and the next level of refinement would then rather 
involve the fluctuations of the vector meson fields. 

The thermal equilibrium values of the two order parameters are then determined by the 
resulting equations of motion for (a) and ((), 



£(°) = Um\ + c — =(^ 2 ) + ^(SK 2 ) + 2(M() (24) 



(± , 

+ (V2\'(5K 2 )-8\(5a5())(() , 



^(C> = V2f K m 2 K - ±=Uml + c ((5n 2 ) + (a) 2 + (6a 2 ) - i(^ 2 >) (25) 
+ (v / 2X'(5K 2 ) -8\(5a8()) (a) , 
where the effective masses for the order parameters are given by 

£ = -/i 2 - 2c(C) + 2A ((5vr 2 ) + (a) 2 + 3(5a 2 ) + l -(5 V 2 ) + (5K 2 )^ , (26) 
+4A' ((5n 2 ) + (a 2 ) + 3(5a 2 ) + (C 2 ) + (5( 2 ) + (5 V 2 ) + (5K 2 )) 

$ = ~f + 4A ((C 2 ) + 3(5C 2 ) + l(Sv 2 ) + <^ 2 >) (27) 
+4A' ((5n 2 ) + (a) 2 + (6a 2 ) + (C) 2 + 3(^ 2 > + (^ 2 ) + (5K 2 )] 



The mixing of the field fluctuations 5a and 5( can be treated to a sufficient degree of 
approximation by the replacements 

(5a 2 ) — > (5a 2 ) cos 2 6 S , 

(5( 2 ) — (5a 2 ) sin 2 9 S , (28) 
(5a5() — > (5a 2 ) sin 6 S cos 9 S , 
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where we have neglected fluctuations of the heavier ( meson. Mixing between 77 and 7/ is 
ignored for the thermal fluctuations as the mixing angle turns out to be near zero even 
at high temperatures. The final expressions for the meson masses can be found in the 
Appendix. 

In order to simplify the notation, we omit the brackets (■) in the following so that the 
order parameters will be denoted by a and ( and the thermal field fluctuations by 5K 2 
and so on. The coupled equations for the order parameter (E4LE51), the meson masses (see 



Appendix), and the fluctuations ( p3| , p8D are solved iteratively with an overdamped Newton 
method. 

We note in passing, that we assume that the coupling constants of the effective interaction 
is temperature independent. This is a strong assumption, especially for the anomaly term 
where one knows that the corresponding coupling constant should be much smaller at the 
chiral phase transition [|50|. Effects should be especially seen for the 77/77' mesons |56| , |57| and 
could even lead to parity- violating effects as pointed out recently |55| . How the coupling 
constants scale with temperature is not precisely known. There exists estimates but they can 
be applied in principle only for a quark-gluon plasma. Therefore, we leave this interesting 
question for future investigations. 



A. Temperature dependence of the order parameters 

Figure |2] shows the temperature dependence of the total fluctuations of the tt, K, a, and 
rj quasi-particle fields for the various values of m a {i.e. 5n 2 = (5ir 2 ) + (Stt%) + (Sn 2 ), etc.). As 
expected, the pion fluctuations are dominant at low temperatures due to their small mass, 
and the other contributions can be safely neglected for T < 50 MeV. The situation changes 
drastically at higher temperatures as the strange degrees of freedom are activated. The 
kaon fluctuations constitute a 20% correction at T = 100 MeV and are already half the pion 
contribution at T = 150 MeV and can thus no longer be neglected. The kaon fluctuations 
become remarkably similar to the ones of the pions at higher temperatures and contribute 
equally at T = 400 MeV. Note that the degeneracy factor for the kaon is four thirds times 
that of the pion so that the kaon contribution would dominate if the masses of the pion and 
the kaon were equal. The thermal fluctuations do not change significantly with m^, except 
of course for those of the a field itself which for large m a grow rapidly in the a crossover 
region as the effective a mass experiences a pronounced minimum (see Fig. ^). 

The figure also displays the thermal evolution of the two order parameters. They exhibit 
a smooth crossover from their large (and nearly equal) vacuum values towards ever smaller 
values at high temperatures. The fastest change for a occurs around T = 210 MeV. The 
decrease of the strange order parameter ( is much less pronounced, as a consequence of the 
larger kaon mass, and the corresponding crossover does not occur until above 400 MeV. 

It is evident that the degree of chiral symmetry restoration achieved depends sensitively 
on the adopted value of m a . For the smaller values the order parameters decrease only 
moderately. For example, at T = 400 MeV a is still above 50 MeV, for the commonly 
employed value of m a = 600 MeV whereas it has dropped below 30 MeV for m a = 1 GeV. 
The more complete symmetry restoration at higher m a values is due to the fact that the 
parameter /z 2 has now become positive and the coupling constant A' has increased. The 
former feature affects the spontaneous symmetry breaking at zero temperature, while the 
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latter enhances the temperature dependence of the meson masses. On the basis of these 
observations, we expect that the most dramatic effects of the extension to SU(3) will occur 
for large values of m a . 

The temperature dependence of the two order parameters is well illustrated by the 
thermal path of the equilibrium point in the (cr, Q) plane, as shown in Fig. [5] for m a = 
600, 800, 1000 MeV. The solid circles mark the progression of the path in temperature in- 
tervals of 100 MeV. In addition, we also show the result of omitting the kaon fluctuations as 
well as the result for the pure SU(2) case. In the latter case, there is obviously no change in 
the strange order parameter and we note that o reaches very small values for temperatures 
above T = 300 (the last two dots are for temperatures of 400 and 500 MeV). The SU(3) 
case without kaon fluctuations is close to the pure SU(2) case but differs from the full SU(3) 
calculation. In particular, ( hardly changes in the absence of kaon fluctuations. The paths 
for the full SU(3) calculation develop very slowly at low temperatures (note that the points 
for T = 100 MeV are quite close to the vacuum point), which reflects the above finding 
that the order parameters change little over that temperature interval. In the intermediate 
temperature range of T = 150 — 300 MeV both order parameters decrease approximately 
linearly and the path follows a nearly straight line. At higher temperatures, a quickly sat- 
urates while ( is only now beginning to drop faster. At these high temperatures, the kaon 
fluctuations now provide the dominant contribution to the meson masses. We see that the 
inclusion of kaon fluctuations (i.e. going from SU(2) to SU(3)) has a dramatic effect on the 
thermal evolution of the order parameters. 

It is interesting to compare our findings for the full SU(3) case with related previous work. 
The linear SU(3) a model was solved in the 1/Nf expansion by means of the saddle-point 



approximation for finite temperatures ||59|| . It was found that the light-quark condensate 
drops sharply at 180-190 MeV, while the strange-quark condensate has changed by only 
20% at T = 300 MeV as compared to its vacuum value. On the other hand, we find that 
the changes of the order parameters depend sensitively on the chosen parameter values, 
especially the value for m a . The heavier the value of m a is, the larger is the drop of both 
order parameters as the temperature is increased and the higher is the critical region in 
which the cross over occurs. The change of the strange-quark condensate shown in Ref. |59| 
is smaller than in our mean-field results, while it is the reverse situation for the light-quark 
condensate. 

There exists only one lattice calculation of the light and strange quark condensates for 
Nf = 2 + 1 [ 53J . The light-quark condensate drops more than 90% through the transition 
region, while the strange-quark condensate decreases only moderately (by 30%) over the 
same temperature region. Our calculations yield a less dramatic behavior, especially for 
a low a vacuum mass, which thus seems to favor a high vacuum value for m a . For the 
case m a = 1 GeV, we obtain a change of only 70% for a and about 30% for (. Still, 
the a condensate changes less than indicated by lattice data, while the strange condensate 
follows the trend seen on the lattice but shifted towards higher temperatures. However, this 
comparison should be taken only as suggestive, as we are using an effective model. 
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B. Meson masses at finite temperature 



The effective meson masses will exhibit a temperature dependence due to the thermal 
decrease of the order parameters and the associated growth of the field fluctuations (the 
explicit expressions for the meson masses, including the contributions from thermal fluctu- 
ations, are given in the Appendix). The effective masses are obtained in connection with 
the above self-consistent solution of the coupled equations carried out for the purpose of 
determining the temperature dependence of the order parameters. 

Figure |] shows the resulting masses for the pseudoscalar nonet (ir, K, r], rf) and the 
scalar nonet (a, a , K* , f ). For temperatures below 80 MeV there is hardly any change in 
the meson masses at all. For intermediate temperatures, 80 < T < 230 MeV, the ir mass 
as well as the K and rj masses start to increase approximately linearly. The rj mass is also 
increasing with temperature but relatively less so than the other pseudoscalar meson masses. 
By contrast, all scalar meson masses decrease with temperature through this range, since the 
decrease of the order parameters bring them closer to their respective critical boundaries. 

For m a = 600 MeV, the a meson mass drops below the K and rj masses and reaches its 
lowest value of 520 MeV at T = 230 MeV. The ao and K* mesons masses change similarly, 
losing about 70 MeV of mass, but the fo mass drops only 35 MeV from its free value. At even 
higher temperatures, T > 230 MeV, the thermal fluctuations dominate all the meson masses 
and they now increase linearly with temperature. The slope for the pseudoscalar mesons is 
steeper than that for the scalar mesons. Especially the pion mass changes remarkably at 
high temperatures, reaching 750 MeV at T = 400 MeV. The a meson becomes degenerate 
with the pion at around that temperature. The K and rj masses are close to one another 
and reach values just below 900 MeV at T = 400 MeV. The rj mass is rising so rapidly with 
temperature that it becomes the heaviest meson of all at temperatures above 400 MeV. 

For m a = 800 MeV, the pseudoscalar meson masses show a similar behavior, except that 
the rj is rising less steeply through the transition region. All scalar meson masses are again 
dropping as the order parameters are dropping towards chiral symmetry restoration, but it 
is now more pronounced. The a mass is falling by nearly 25% to m a = 615 MeV at T = 
270 MeV. The masses of ao and K* drop from their free values by 100 MeV. Furthermore, 
the fo loses about 85 MeV, but the minimum is reached at the higher temperature of 
T = 370 MeV. Again, the ir and a masses, as well as those of rj and K, are all becoming 
very similar behavior above 300 MeV. 

The changes for the scalar meson masses are even more pronounced for m a = 1 GeV. 
The a mass now drops by 28% and reaches its minimum value of 716 MeV at T = 305 MeV. 
The a and K* are dropping 125 MeV and 146 MeV, respectively. The minimum of the K* 
mass is located at a slightly higher temperature, T = 325 MeV, due to its strong coupling 
to the strange order parameter ( which starts decreasing at a higher temperature than a. 
The fo now shows a significant mass change, having decreased from 1550 MeV to 1240 MeV 
at T = 400 MeV. If continuing to higher temperatures, a minimum mass of 1220 MeV 
would be reached at T = 440 MeV. Here the fo mass is strongly influenced by the strange 
order parameter, hence the mass minimum is shifted to higher temperatures as compared 
to the other scalar mesons. The pseudoscalar meson masses are showing a steeper rise with 
temperature, as compared to the other cases with a lower vacuum value of m a . The rj mass 
shows a different behavior, as it remains rather unchanged up to T = 300 MeV, before 
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starting its growth. At T = 400 MeV, the masses of K, 77, ao, and K* are closely clustered 
around 1 GeV. The rf and fo masses appear to converge at higher temperatures, signaling 
the onset of full chiral symmetry restoration. 

Comparing the results obtained with the different a masses, we observe some general 
trends. First, the pions, as well as the other pseudoscalar mesons, grow more massive at 
high temperatures; there are thus fewer medium-modified pions in a thermalized system 
at high temperatures as compared to a free gas at the same T. Second, the effective K 
and 77 mesons become approximately degenerate with the effective 7r meson. This will 
enhance the yield ratios to K/tt ~ | and r]/n « | as compared with the free-gas estimate 
of K/n « |e" Am / T w 0.5 and r]/ir « 0.1 for T = 400 MeV. 

At high temperature, the thermal fluctuations of the vector mesons can no longer be 
neglected. Of course, their impact on the results presented depend on the change of their 
masses with temperature. If the vector meson masses stay near their vacuum values up to 
T = 300 MeV, or perhaps drop with temperature, the corresponding thermal fluctuations 
will grow important already in the transition region for the order parameter a. On the 
other hand, if the vector meson masses increase with temperature, as does the 77', then their 
contributions can possibly be neglected over the temperature range shown, as they are too 
heavy to become significantly agitated. 

The thermal ao and K* fluctuations become important at high temperatures, say T > 300 
MeV, as their masses are then close to those of rj and K, respectively. As the thermal masses 
are in the range of 1 GeV, vector mesons should be also included and might constitute a 
bigger correction term at very high temperatures. Therefore, we omitted the ao and K* 
fluctuations. It might be interesting to study vector mesons within the linear a model and 
they should clearly be included in studies of the high-temperature behavior of pseudoscalar 
and scalar mesons above the chiral transition. 

A more important effect is associated with the temperature dependence of the anomaly 
term (Uf). If the anomaly term becomes considerably weaker as the temperature is increased, 
it will reduce the thermal masses of the rj' rather than increase them. The thermal mass of 
the a mesons will also be reduced by this effect. First estimates indeed show that this is the 
case, so that these thermal fluctuations then ought to be included. Whether the presence of 
these additional low-lying mass states will change the order of the phase transition is beyond 
the scope of the present investigation but will be addressed in a forthcoming work. 

V. DYNAMICS 

In this section, we consider a simplified dynamical scenario that approximates the con- 
ditions occurring in the central region of a relativistic heavy-ion collision. For this purpose, 
we adopt the approach given in Ref . |22| , in which a scaling expansion is emulated by means 
of a Rayleigh cooling, and extend it to SU(3). We discuss the time evolution of the order 
parameters as well as that of the pseudoscalar and scalar meson masses. 
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A. Scaling expansion in SU(3) 



As an initial state, we consider thermalized matter at a specified temperature, with only 
the scalar and pseudoscalar mesons included. The matter is then assumed to experience 
a uniform scaling expansion which causes the field fluctuations to decrease according to 
the dimensionality D of the expansion. This brings the system out of equilibrium, as the 
effective potential governing the order parameter reverts to its vacuum form faster than 
the order parameter can follow. For the standard SU(2) case, a boost invariant expansion, 
corresponding to D — 1, was considered in Ref. while the three-dimensional scenario 
was explored in Refs. jn^jlj. If the expanding system is strictly uniform, its properties 
depend only on the elapsed proper time r, i.e. the time experienced in a local comoving 
frame, and the dAlembert operator is then given by 

J_d_ D d_ <P_ Dd_ 

t d dr dr dr 2 t dr 

The equations of motion of the two order parameters are therefore readily obtained, 
/ + — + = Uml + c ( -^5rf + ^5K 2 + 25a5(\ + 



\dr 2 rdr ™T J *""* "\3V2 ' ^2 s ) 

[V2X'5K 2 - 8\(5a5()) ( (30) 

+ ~ + ^ C = V2f K m% - -±=Uml + c (V + a 2 + 5a 2 - i V) 

+ (V2X'5K 2 - 8\8a8() a , (31) 



where the effective masses and /i^ are given in Eqs. ( p6|) and (|27|). 



The equations of motion (|30| - |3T1) , can be used to propagate the order parameter. At 
each time step, the effective masses are obtained by means of the expressions given in the 
Appendix. In order to approximately emulate the effect of the expansion, we assume that 
the field fluctuations exhibit a simple power behavior, 

5n(r) 2 = 5n(r r , 5K{rf = 8K{r f {^f , 

^(r) 2 = ^(r ) 2 (^) D , 8*(rf = 6a(r f , (32) 

as they will eventually when the system has become so cold that the individual quasi-particle 
modes are decoupled (to denotes the initial time when the expansion is started). Because of 
the mixing between 5a and 5(, and the interdependence of the various quantities, a simple 
iteration procedure is carried out until convergence is achieved at each time step. 



B. Evolution of the order parameter 

In our first illustrations, we employ a relatively high initial temperature, T = 400 MeV, 
in order to achieve a high degree of chiral symmetry restoration in the initial state. Moreover, 
the initial time is chosen as tq = 1 fm, as is most commonly done, and &{tq) = and 



14 



C( T o) — 0- Figure [5] shows the trajectories of the two order parameters for m a = 1 GeV 
and in two different expansion scenarios, a longitudinal (D = 1) and isotropic (D = 3). 
The former represents a lower bound on the degree of expansion, as it assumes that no 
transverse expansion occurs at all, whereas the latter probably overestimates the expansion, 
so the actual physical scenarios is expected to be intermediate between these two extremes. 

For the three-dimensional expansion, the strange order parameter is changing first and 
overshooting its vacuum value. The a order parameter remains quite unchanged during that 
time and starts growing towards its vacuum value only at later times. The long-dashed curve 
traces the movement of the equilibrium point as the temperature is reduced; this thermal 
path would be followed if the expansion were adiabatic so the system had time to equilibrate 
throughout the expansion process. However, because the expansion is fast, the dynamical 
trajectory oscillates around the thermal path. For the fastest expansion, only one oscillation 
occurs before the effective potential has reestablished its vacuum form and the remaining 
oscillations are then around the vacuum point. It is interesting that the early oscillations of 
( are much stronger than those of a, whereas the situation is reversed for the late oscillations 
around the vacuum point. 

The one-dimensional expansion shows a distinctly different behavior. As above, the 
initial evolution is primarily in the direction of the strange order parameter, but there are 
now several oscillations around the adiabatic path before the vacuum point is reached. This 
is because the order parameter trajectory oscillates around the instantaneous minimum in 
the effective potential which now develops relatively slowly from its initial value near the 
origin towards the vacuum point. Thus there is more opportunity for the order parameter 
to overshoot the vacuum point and rebound towards smaller values. 

C. Evolution of the thermal fluctuations 

The pattern in the evolution of the order parameters will be reflected in the behavior of 
the pseudoscalar and scalar mesons masses. Figure |j] (left) shows the time dependence of the 
pseudoscalar meson masses for the one-dimensional expansion. The masses are scaled to their 
respective vacuum values so that they will tend to unity at large times. The pseudoscalar 
meson masses remain well above zero during the evolution. The initial value of the pion 
mass is much larger than its vacuum value and it exhibits a slowly decreasing average 
behavior, overlayed with large oscillations containing several different frequencies. This 
rather complicated pattern is a result of the coupled evolution of the two order parameters 
during the expansion. The oscillations are nearly undamped over the entire time interval 
shown, i.e. up to r = 15 fm. The K and rj masses also show oscillations but they appear 
to have only one main frequency and the amplitudes are much smaller than for the 7r mass. 
The r)' mass does not change very much and has only small fluctuations around its vacuum 
value. 

The scalar masses experience a different time evolution, as depicted in Fig. |] (right). 
All scalar masses start below their vacuum value and approach their free value from below, 
contrary to the behavior of the pseudoscalar mesons. The oscillations of the scalar meson 
masses are sometimes in and sometimes out of phase. Especially the a and the ( meson 
masses are quite often in opposite phases. The amplitudes of the oscillations are rather 



15 



similar, with that of the a meson mass being the largest, and they are not sufficiently large 
to overshoot the respective vacuum value. 

Performing a Fourier transform of m 2 (r), we find that there are two main frequencies 
involved, one at 0.74 fm -1 and the other at 1.2 fm _1 . The oscillations of the a order 
parameter is dominated by the lower frequency, with a weaker component from the higher 
one, whereas the strange order parameter has only the higher frequency. The pion has both 
frequencies with the lower one being stronger. The K, rj, and K* masses follow the frequency 
of the strange order parameter, while those of rf and seem to have both frequencies at 
about equal strength. 

Now we turn to the three-dimensional expansion. Figure [7] (left) shows the time depen- 
dence of the pseudoscalar masses. The oscillations of the pion mass are now much stronger 
than for the other pseudoscalar mesons (and the pion curve is therefore scaled down by a 
factor of ten). The oscillations are also much smoother than in the one- dimensional case. 
Most significantly, the oscillations brings the pion mass into the region of instability sev- 
eral times (a total of six in the specific calculation). Remarkably, the second minimum 
at t pd 3 fm is slightly deeper than the first minimum. This is a result of the dynamical 
interplay between the two order parameters. Indeed, the order parameters are not really in 
phase at the time of the first minimum, as the initial change of ( is much larger than that 
of a. However, at the time of the second minimum the two order parameters are already 
oscillating around the vacuum point and are in phase, so that the oscillations in the pion 
mass are enhanced. At r « 3 fm/c the masses of all the other pseudoscalar mesons also 
exhibit their lowest minimum, but the mass oscillations of these mesons are weaker than in 
the one-dimensional expansion (note the different mass scales in Figs. |6| and [7]). The K and 
r) masses are evolving very similarly, but the rf mass shows a different behavior: it is not 
oscillating with one dominant frequency, as was the case for the tt, K, and rj masses. For 
all pseudoscalar mesons, the damping of the oscillations is stronger than in the longitudinal 
expansion, as would be expected. 

The scalar meson masses obtained for D = 3 are plotted in Fig. [Tj (right) and their 
oscillations are seen to be stronger than for D — 1. Again, the scalar meson masses start 
below their vacuum values but then overshoot them later. The a meson mass has a deep 
minimum at r = 1.4 fm/c which comes remarkable close to zero. Also the ao has strong 
oscillations in phase with the pion. The frequencies of the K and ( masses differ from those 
of the it and rj masses. 

A Fourier decomposition of m 2 (r) reveals that there are again two main frequencies, 
0.78 c/fm and 1.22 c/fm. The tt, a, and ao have the lower frequency and the K* and ( have 
the higher one, while the rf has both frequencies with equal strength. By contrast, the K 
and rj masses evolve with a broad distribution of frequencies around 1.3 — 1.56 c/fm. 

Finally, we note that since the displayed masses are scaled relative to their free values, the 
pion will always seem to have the strongest oscillations which could be somewhat misleading. 
In fact the pseudoscalar mass shifts are all of the same order of magnitude, around 150- 
400 MeV for the first extremum in the evolution for D — 3. On an absolute scale, the 
fluctuations of especially the K and 77 masses are as large as those of the tt mass which my 
possibly lead to observable effects in relativistic heavy-ion collisions. 
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D. Spectral enhancements 



We finally utilize the calculated time-dependent effective ir and K masses to illustrate 
the possible associated enhancements in the corresponding spectral distributions of the free 
mesons. Towards this end, we assume that the evolution of the field amplitude for a given 
meson mode is governed by a Klein-Gordon field equation, 

0k = u 2 k k , u 2 k = m 2 (t) + k 2 . (33) 

Here k denotes the momentum of the mode and the time-dependent frequency is given by the 
indicated dispersion relation, where m 2 (t) is the calculated time-dependent effective mass 
squared. Invoking the method presented in Ref . |H] , we can then calculate the enhancement 



coefficient for that particular mode, X^. This quantity (which is never less than unity) 
expresses the degree to which the population of the mode has been enhanced due to the 
time dependence experienced by the frequency. Thus, if the initial occupation number is 
z/k(ro) = ( 7 "o| a k a k| r o)) then its final value is given in terms of X^ as 

M r f) = ( T f\ a W\ T f) = Mr ) + |]X k -|, (34) 

where the presence of the one half reflects the effect of the quantum fluctuations of the field 
amplitude for that mode. 

Figure |8] shows the resulting pion enhancement coefficients for the two pseudo-expansion 
scenarios considered. For D — 1 there is a significant enhancement of the softest modes 
as well as a moderate peak near a pion momentum of 400 MeV. The former feature is a 
direct reflection of the relatively slow overall decrease of m 2 in the course of time, while 
the peak near a momentum of 400 MeV/c results from the quasi-periodic oscillations of the 
frequency. This latter effect was discussed already in connection with the enhancement of 
electromagnetic processes The results for D = 3 exhibit a large and fairly sharp peak 

at a slightly higher momentum, with a more moderate and smooth enhancement of the soft 
modes. The latter reflects the fact that m 2 (t) remains fairly constant on the average, after 
the initial fast drop. As a consequence of the more perfect quench achieved with D = 3, the 
oscillations of m 2 are more harmonic and, therefore, the parametric amplification becomes 
stronger. The resulting peak in the enhancement is located at a momentum of 480 MeV 
which corresponds very closely to the expected resonance energy, u)\ = \m a = 500 MeV. 

Figure ^| shows the corresponding results for the kaons. Since the relative modification 
of the effective kaon mass is smaller (see Figs. ||]7]), the associated enhancements are also 
much more modest. Nevertheless, it is possible to discern features similar to those of the 
pions. In particular, for D = 1 there is relatively broad enhancement of the softer modes 
followed by a small resonance peak and for D = 3 the most prominent feature is resonance 
at a slightly larger momentum but, perhaps surprisingly, it is quite small. 

We finally wish to stress that these illustrations are based on a fairly schematic picture 
and they should not be taken as quantitative predictions. Nevertheless, they do suggest 
that significant spectral enhancements may be present in the SU(3) description and they 
also give a general idea of their specific features. Evidently, a quantitative prediction would 
require much more elaborate calculations of the complicated collision dynamics. 
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VI. CONCLUSIONS 



The present investigation presents a first exploration of how the inclusion of strangeness 
in the effective a model affects the features associated with disoriented chiral condensates. 
The resulting SU(3) model contains both the scalar and pseudoscalar meson nonets and 
introduces the strange order parameter, ( = (ss), thus making both the phase structure and 
the dynamics more complex. 

The inclusion of strangeness is expected to have significant effects once the temperature 
exceeds the mass of the strange quark and thus it should not be ignored in DCC discussions. 
Indeed, the crossover towards chiral symmetry is being pushed to higher temperatures, even 
for the non-strange order parameter a, though the equilibrium value of o is reduced when 
the a mass is increased. Consistent with this trend, the thermal masses of the 18 mesons 
then also become more similar. 

By considering suitable idealized expansion scenarios, we have sought to ascertain the 
effect of the model extension on the non-equilibrium dynamics following a high-energy nu- 
clear collision. The field fluctuations reflecting the presence of kaons are important for the 
dynamics and can not be neglected. Furthermore, these calculations yield significant spec- 
tral enhancements not only for the pions (as was the case in the SU(2) description) but also 
for the kaons (though to a much smaller degree). 

We have pointed out that the neutral fraction R K of kaon emitted by a single DCC 
domain has a constant distribution, Px(Rk) — 1 which is qualitatively as anomalous as 
the inverse square-root distribution P 7T (R n ) = Xjl^fR^ governing the neutral fraction of 
DCC pions. Relative to this latter distribution, the kaon distribution P K (R K ) may be more 
amenable to measurement by suitable event-by-event analysis due to the fact that Kg decays 
into charged pions. 

Finally, if chiral symmetry is approximately restored, the K/n and f]/ii equilibrium 
ratios are enhanced by factors of 2-3 due to the mass shifts of the mesons. If the expansion 
is sufficiently fast, this effect may be observable, even though the signal will be reduced by 
decaying resonances. 

The present investigation assumes that the coupling constants are temperature inde- 
pendent. A possible temperature dependence of the anomaly term will certainly change 
the meson mass spectrum at the higher temperatures, as well as the evolution of the order 
parameters and the effective masses. This will be studied in a forthcoming work. 

Our present exploratory study suggest that various effects associated with the chiral 
phase transition in SU(3) might be observable at RHIC and we hope that our present 
findings will stimulate more elaborate investigations. 
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APPENDIX 

The effective masses for the pseudoscalar mesons ir, K, rj, and a are given by 

ml = -fx 2 + 2A (Jtfvr 2 + SK 2 + Sri 2 + a 2 + Sa 2 

+ 4A' (Jc^tt 2 + 5K 2 + 5 V 2 + a 2 + bo 2 + C 2 + 5C 2 

-li 2 + 2A ^tt 2 + hiX 2 + 5rf + a 2 + 5a 2 + 2( 2 + 25( 2 - V2 {a( + 5a5C)j 

+ 4A' ^vr 2 + ^-bK 2 + 5r, 2 + a 2 + Sa 2 + ( 2 + 6( 2 
— V2ca , 



m\ 



m P0Q 



V + 4A (V + 6K 2 + 5r 1 2 + l - (a 2 + 5a 2 + C + K 
+ 4A' (Stt 2 + SK 2 + Sr, 2 + a 2 + Sa 2 + ( 2 + 6( 2 ) 
+ (V2a + C) , 

-ll 2 + 2A Un 2 + 6K 2 + 35r] 2 + - (a 2 + 5a 2 + 4( 2 + 45( 



3 

4A' (6n 2 + 5K 2 + 36 V 2 + a 2 + 5a 2 + C 2 + S( 2 
2 



m P88 



--c(2V2a-() 

nf lw = V2\ (25n 2 - 5K 2 - 25r] 2 + | (a 2 + <5a 2 - 2( 2 - 25( 2 

-L(a-V2C) , (35) 



where the thermal fluctuation of the pion field is denoted by Sw 2 and similarly for the other 
fields. The mixing between the thermal fluctuations of rj and rj' are ignored because of their 
small mixing angle. For finite temperature, one gets also a mixing between ir and rj due to 
the thermal kaon fluctuations. This term vanishes in the mean-field approximation in the 
above expressions but it would give a mixing term ml vo = 2y/3X5 K 2 . Since this term will 
be important only at high temperatures, where the 7r and rf masses are closer in value and 
the thermal kaon field fluctuations are sizable, we ignore it here. 

Using the notation of the Particle Data Group |39[], we define the pseudoscalar mixing 
angle as 



rf \ _ ( cos 6 P sin 6 P \ ( r) 
r] ) V - sin p cos 6 P J \t] 8 



(36) 



The physical masses of the rj and rf are given by the mixing formula j62[, 



19 



m »7 = m2 P88 cos2 @p + m Poo s i n2 @p ~~ m P08 sm , (37) 
m 2 , = mp 88 sin 2 6* p + m 2 P00 cos 2 # p + mp 08 sin 2# p , (38) 



where the pseudoscalar mixing angle 9 P is fixed by 

tan2fl p = 2 2m ™ 2 • (39) 

m P00 m P88 

These mixing formulas relate in a unique way the eigenvalues of the mixed mass matrix to 
the physical unmixed masses m v and m v i. Note that the physical i] state is then the one 
closest to the octet state rjs, while the 7/ is the state closest to the singlet state 770- When 
the masses are defined in the above way, then n and rj' are degenerate when the anomaly 
term vanishes (c = 0). 

The corresponding expressions for the scalar mesons are 



ao 



V + 2A (W 2 + 5K 2 + Urf + 3a 2 + 35a 2 ^j 
+ AX' (5tt 2 + 5K 2 + 5 V 2 + a 2 + 5a 2 + ( 2 + 5( 2 ) 

-fi 2 + 2A (5n 2 + 5K 2 + 7 -5rf + a 2 + 5a 2 + 2( 2 + 25( 2 + V2(a( + 5a5()j 

+ AX' (5n 2 + 5K 2 + 5 V 2 + a 2 + 5a 2 + ( 2 + 5( 2 ) 
+ V2ca , 



m 500 — 



-/i 2 + 2A (5-k 2 + 5K 2 + ^5 V 2 + 3 (a 2 + 5a 2 
+ AX' (5n 2 + 5K 2 + 5rf + 3a 2 + 35a 2 + C 2 + S( 2 ) 



2< , 

2 



m S88 



m S08 



-fi 2 + 4A (W 2 + ^ 2 + 3C 2 + 35C 2 ) 
+ 4A' (5tt 2 + 5K 2 + 5 V 2 + a 2 + 5a 2 + 3( 2 + 35( 2 ) , 

-V2X5K 2 + 8 A' (crC + 5a5() - 2ca . (40) 



The physical masses for the scalar mesons a and ( are obtained in analogy with the 77 and 
7/ described above, where 



a \ _ ( cos 6 S sin 6 S \ f a 
CJ \-sin6> s cos9 s ) \<j 8 



(41) 



with 9 S being the scalar mixing angle. 
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FIGURES 
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Order parameter a [MeV] 



FIG. 1. The critical boundaries for the pseudoscalar mesons in the plane of the order parameters 
a = (qq) and £ = (ss). The cross indicates the vacuum equilibrium point (oVaoCvac) wn ere all 
meson masses are positive. The results were obtained with four different values of m CT as indicated. 
The curves delineate the various critical boundaries, where the respective effective mass vanishes. 
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FIG. 2. The temperature dependence of the it, K, rj, and a field fluctuations in thermal equi- 
librium, together with the equilibrium values of the two order parameters a and £, for m a = 400, 
600, 800, and 1000 MeV. 
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FIG. 3. The thermal path of the two order parameters for various cases of interest, as indicated. 
The solid dots are plotted in temperature steps of 100 MeV. All paths start from the vacuum point 
at the upper right corner and the points for T = 100 MeV are still very close to the vacuum point. 
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FIG. 4. The temperature dependence of the effective masses of the pseudoscalar (-7T, if, 77, 7/) 
and scalar (<r, if*, ao) meson nonets for the four values of m a studied, as indicated. 
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110 




FIG. 5. The dynamical path of the order parameter (cr, £) as a result of a uniform 
pseudo-expansion in one and three dimensions, starting from thermal equilibrium at T = 400 MeV. 
Time steps with At = 0.2 fm are indicated along the trajectories with squares (D=l) and circles 
(D=3). The thermal path is shown for comparison (long-dashed curve) and the cross indicates the 
location of the vacuum, (a vac , Cvac), towards which all trajectories converge in time. 
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FIG. 6. The time evolution of the (squares of the) effective masses (a) for the pseudoscalar 
mesons n, K, rj, rj' and (b) for the scalar mesons a, £, K* , ao in a pseudo-expansion with D = 1. 
The masses are given relative to their respective free values and m a = 1 GeV has been used. 




FIG. 7. The time evolution of the (squares of the) effective masses (a) for the pseudoscalar 
mesons ir, K, rj, rf and (b) for the scalar mesons a, (, K* , ao in a pseudo-expansion with D = 3. 
The masses are given relative to their respective free values and m a = 1 GeV has been used. 
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FIG. 8. The spectral enhancement coefficient for pions as obtained on the basis of the 
calculated time-dependent effective masses for the two pseudo-expansion scenarios considered, in 
both cases starting from matter in thermal equilibrium at T = 400 MeV. 
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FIG. 9. The spectral enhancement coefficient for kaons, as obtained on the basis of the 
calculated time-dependent effective masses for the two pseudo-expansion scenarios considered, in 
both cases starting from matter in thermal equilibrium at T = 400 MeV. 
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TABLES 



TABLE I. The meson masses in vacuum for the three parameter sets corresponding to 
m a = 400, 600, 800, and 1000 MeV. The pseudoscalar masses, the pseudoscalar mixing angle, 
and the masses of ao and K* are independent of that choice and are given in the first line. The 
coupling constants and the masses of the a and £ mesons with their mixing angle 9 S are listed. All 
masses are given in MeV. 
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